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A GENERAL FORMULA FOR THE ALGEBRAIC DEGREE IN 
SEMIDEFINITE PROGRAMMING 

HANS-CHRISTIAN GRAF VON BOTHMER AND KRISTIAN RANESTAD 

Abstract. In this note, we use a natural desingularization of the conormal variety of the 
variety of (n x n)-symmetric matrices of rank at most r to find a general formula for the 
algebraic degree in semidefinite programming. 



1. The algebraic degree in semidefinite programming 

Let P m be a general projective space of symmetric (nxn)— matrices up to scalar multiples, 
and let Y r C P m be the subvariety of matrices of rank at most r. In this note we find a 
formula for the degree 5(m, n, r) of the dual variety Y* whenever this is a hypersurface, i.e. 
5(m, n, r) is the number of hyperplanes in P m in a general pencil that are tangent to Y r at 
some smooth point. 

In [3] the algebraic degree of semi-definite programming is introduced and is shown to 
coincide with a degree of a dual variety as above j3l Theorem 5.2]. In particular, the 
number 8(m, n, r) is a measure of the algebraic complexity of the rank r solution for a 
general objective function optimized over a generic m-dimensional affine space of symmetric 
(n x n)— matrices. 

Our general formula for S(m, n, r) extends the results of [3J. It is expressed in terms of a 
function on subsequences of {1, n}. Let 

= 2* -1 , if) it j = ^ [ + ^ 2 ) when * < 3> 

k=i ^ ' 

and 

Ai,...,i r = P£(*Pi k ,ii)i<k<i<r if r is even 
■0ii,...,v = Pi(ipi k ,n)o<k<i<r if r is odd 
where ipi 0) i k = ipi k and Pf denotes the Pfaffian. 

Theorem 1.1. The algebraic degree 

5(m,n,r) = y^i^jc, 
i 

where the sum runs over all strictly increasing subsequences I = ...,z n _ r } of {1, ...,n} 
of length n — r and sum %\ + ... + i n _ r = m, and I c is the complement {1, n} \ I. 

Notice that the formula makes sense only when m > ( n ~l +1 ) and ( r ^ X ) < ( n ^ 1 ) — m, 
which is precisely the allowable range for the algebraic degree (cf. [3], Proposition 3.1]). 
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2. The bidegree of a conormal variety 

Let V and V* be an n- dimensional vector space and its dual with fixed dual bases. Then 
we interpret a symmetric (n x ri)— matrix as a bilinear map <$$ '■ V —>■ V* (or (fix '■ V* — > V) 
such that (j)(vi)(v2) = 0(1*2) (^1) for any pair of vectors t>i, t>2. The symmetric maps V — > V* 
form S" 2 ^*, while the symmetric maps V* — > V form S' 2 V. In this notation we let F r C 
P(S' 2 V A ) be the space of symmetric (n x n)— matrices of rank at most r (up to scalars) and 
CY r C F(S 2 V*) XqF(S 2 V) the variety of pairs of symmetric matrices (S, T) of rank at most 
n — r and r whose matrix product ST = 0. The variety CY r is the conormal variety of Y r , 
i.e. the closure of the incidence set of pairs (S,T) such that S G F(S 2 V*) corresponds to a 
hyperplane in F(S 2 V) that is tangent to Y r at the smooth point T. 

The bidegree of CY T C P(S' 2 l A *) Xc P(S' 2 V r ) is the generating function for the algebraic 
degree. More precisely, denote by H\ be the class of a hyperplane on F(S 2 V*) and Hi be 
the class of a hyperplane on P(S' 2 l / ), then the class of CY r is a homogeneous polynomial in 
Hi and i?2- 

Proposition 2.1. [3J Theorem 4.2] 

[C%] = ^i( m ,n,r)FP" m F 2 m 

m 

We may therefore compute the algebraic degree on CY r : 
Corollary 2.2. 

5{m, n,r)= [ H™~ 1 Hp^' 1 n CY r . 

JCYr 

3. A DESINGULARIZATION OF THE CONORMAL VARIETY 

Following the argument of [3J Lemma 3.2] we define a natural desingularization of CY r : 
Fix a rank r subspace K of V. The set of pairs (S, T) of symmetric bilinear maps S : V V* 
and T : V* — > V such that K C kerS* and K 1 - C kerT C V* form a product of projective 
spaces: S determines an element in S 2 (V/K)* C S ,2 V*, while T determines an element in 
S 2 K C S 2 V, since (K*/^) is naturally isomorphic to K*. 

Varying K we find a desingularization of CY r as a fiber bundle over the Grassmannian 
of rank r subspaces of V. More precisely, consider the Grassmannian G = G(r, V) of rank 
r subspaces, and the universal sequence 

on G. Consider the symmetric squares S 2 U and S 2 Q* and their projectivisations ¥(S 2 U) 
and P(S 2 Q*). The inclusions S 2 f/ C S 2 V and S 2 Q* C SV* define a natural map 

X := F{S 2 Q*) x G F{S 2 U) -> P(5 2 ^*) x c P(5 2 K). 

Proposition 3.1. let G = G(r,y) and X = F(S 2 Q*) x G F(S 2 U), where Q and U are 
the universal quotient and subbundle on G. Then the image of the natural projection X — > 
F(S 2 V*) x c P(S' 2 V A ) the conormal variety CY r of the variety Y r of symmetric bilinear 
maps T : V* —>■ V of rank at most r, and X — > CY r is a desingularization of CY r . 
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Proof. On the one hand, the image of the projection consists of all pairs (S,T), such that 
the composition ST = and the ranks of S and T are at most n — r and r, so the image 
is precisely the co normal variety CY r of Y r . On the other hand the variety X is clearly 
smooth, and the map X — > CY r is birational, so it is a desingularization of CY r . □ 

4. The proof of the formula 

We denote by Hi and H 2 also their pullback to X. It remains by Corollary [2^2] to compute 
EIR\ on X, when a + b = - 2= r(n-r) + f^ 1 ) - 1 + ( n ' r 2 +1 ) - 1 = dimX. Let 

Pi : F(S 2 Q*) — > G and p 2 : F(S 2 U) — ► G be the natural projections and p : X — > G be the 
induced projection from the fiber product. We use Segre classes to move the computation 
from X to G. Our projective bundles are the varieties of rank one subbundles and not rank 
one quotient bundles as in [T], so the Segre classes in [fl 2.4] are the Segre classes of the 
dual bundles in our notation. Since S 2 Q* has rank ( n ~l +1 ) , we therefore get 

(pO # (#f n x) = (-iy Si s 2 Q* n G = Sl s 2 Q n G. 

Thus 

(piUH? nx) = s a _ {n -r +1)+1 s 2 Q n G. 

Similarly we compute iff: 

( P2 )M n = s b ^ r)+1 s 2 u* n G. 

On the fiber product we get 

p*(^ a ff 2 6 ni)= p^^ifl n = s a _^-r +1)+1 (s 2 Q)s b _ { r tl)+1 (s 2 u*) n G. 

Combined with Corollary 12.21 we have therefore reduced the computation of S(m,n,r) to a 
calculation on the Grassmannian variety G. 

Proposition 4.1. The algebraic degree 

5(m,n,r)= / s m _^~r+i^(S 2 Q)s^+iy m _^+i^(S 2 U*) f] G. 

Proof. Since X is birational to CY r and H™ H 2 2 fl GF r is finite, 

5(m, n, r) = / H^H^ n GF r = / H^H^^ 1 n X. 

□ 

In this formulation, computing the degree involves Schubert calculus on G. In two steps 
we will circumvent this, following some basic results on Schur functions on Grassmannians 
(cf. |2J and [lj) : First, we express the total Segre classes of the symmetric product with 
the Segre classes Sj = Si(E) of the original bundle: 

s(S 2 E) = Y,^iSi(E) 
i 
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where / runs over all strictly increasing sequences ...,i e } of positive integers of length 
equal to the rank e of E, and sj(E) is the Schur function 

/ s^-x s i2 _i ... S ie -i\ 



si(E) = det 



•Sii-2 s i 2 -2 



Sie-2 
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e ^12— e ■■■ u ie—e y 

In particular sj(E) has degree %\ + ... + i e — ( e 1 2 1 )- The coefficient ipj is the integral valued 
function 



where J runs over all strictly increasing sequences of length equal to the length of J, and 
Ai j is the minor of following matrix with rows indexed by I and columns indexed by J 



A 



We refer to |TJ, Proposition A. 15], for the following closed formulas for the function ipj: 
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and 



where ip. 
Thus, 



fa, 



E 

k=i 



i + 3 ■ 
k 



,ir ~ Pf (V'jfc,ii)l<fc</<r 
,ir = ('0ifc,ii)o<fe</<r 



«0>«fe 



ipi k and Pf denotes the Pfaffian. 



E 

i 



when i < j, 

if r is even 
if r is odd 



4>isi(Q) 



where / runs over all strictly increasing sequences of positive integers / = {ii,i2, ■■■■,in-r} 
with sum i\ + ... + i n _ r = m. (Notice that our sequences are shifted once to the right 
compared to those of [I], in particular they do not contain 0.) 
Secondly, if / has length n — r and J has length r, then 

sj(U) Sl (Q) = SJI V G 

where JI is the concatenated sequence (cf. [H 3.1]). Now, s ji = whenever the sequence 
JI has repeated entries. Furthermore, since Vq is a trivial bundle, sjiVg fl G = whenever 
the degree of sjiVq is positive, i.e. the sum of the entries in JI exceeds ("^ x ) . We conclude 



that sjiVoHG is nonzero only if JI is a permutation of {1, 
and, by (cf. [1, 3.2.1]) 

8i(Q)sj(U*) nG = 5 /JC 



, n} . In particular, i„_ r , > < n 
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where / and J are subsequences of {1, ...,n} of length n — r and r respectively, and the 
sequence J c is the complement {1, ...,n} \ J and, finally, 5ab = 1 if A c = B and 5ab = 
otherwise. 

We put the two steps together and the Theorem 11.11 immediately follows. 
Example 4.2. We compute 5(m,n,r), when n = 5. The relevant values of ip are: 

7 1 2 3 4 5 
ipr 1 2 4 8 16 

/ 1,2 1,3 1,4 1,5 2,3 2,4 2,5 3,4 3,5 4,5 
1 3 7 15 3 10 25 10 35 35 



/ 1,2,3 1,2,4 1,2,5 1,3,4 1,3,5 1,4,5 2,3,4 2,3,5 2,4,5 3,4,5 
V>j 1 4 11 6 23 27 4 18 30 20 

/ 1,2,3,4 1,2,3,5 1,2,4,5 1,3,4,5 2,3,4,5 
t/ji 1 5 10 10 5 

So 

5(1, 5, 4) = ^1^2,3,4,5 = 5 5(2, 5, 4) = ^1,3,4,5 = 20 
5(3, 5, 4) = ^1,2,4,5 = 40 5(4, 5, 4) = ^1,3,4,5 = 40 
5(5, 5, 4) = ^1,2,3,4 = 16 5(3, 5, 3) = ^3,4,5 = 20 

5(4, 5, 3) =^1,3^2,4,5 = 90 

5(5, 5, 3) = ^1,4^2,3,5 + ^2,3^1,4,5 = 7 ■ 18 + 3 ■ 27 = 207 
5(6, 5, 3) = ^1,6^2,3,4 + ^2,4^1,3,5 = 15 ■ 4 + 10 ■ 23 = 290 
5(7, 5, 3) = ^2,5^1,3,4 + ^3,4^1,2,5 = 25 ■ 6 + 10 ■ 11 = 260 
5(8, 5, 3) = ^3,5^1,2,4 = 140 5(9, 5, 3) = ^,5^1,2,3 = 35 
Example 4.3. With a computer it takes less then one minute to compute 

5(105,20,10) = 167223927145503062075691969268936976274880 
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